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To extend the applicability of density functional theory for superconductors (SCDFT) to systems 
with significant particle-hole asymmetry, we construct a new exchange-correlation kernel entering 
the gap equation. We show that the kernel is numerically stable and does not diverge even in the low 
temperature limit. Solving the gap equation for model systems with the present kernel analytically 
and numerically, we find that the asymmetric component of electronic density of states, which has 
not been considered with the previous kernel, systematically decreases transition temperature (T c ). 
We present a case where the decrease of T c amounts to several tens of percent. 

PACS numbers: 74.20.-z, 74.62.Dh, 74.62.Yb 



I. INTRODUCTION 

Non-empirical calculation of superconducting transi- 
tion temperature (T c ) has been a great challenge in theo- 
retical physics. Conventionally, T c s of phonon-mcdiatcd 
superconductors arc estimated by the McMillan- Allen- 
Dynes formula^ 
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where ui\ n , A and /i* (Ref. Q) are the logarithmically 
averaged phonon frequency, electron-phonon coupling, 
and Coulomb pseudo-potential, respectively. This for- 
mula clearly indicates the extreme sensitivity of T c to 
the phonon and electron properties; to achieve a reliable 
calculation of T c , we need to evaluate these factors with 
high accuracy. For the phonon properties, developments 
in methods based on density functional theory^ enables 
one to obtain reliable values very efficiently. On the other 
hand, for while some calculation based on the static 
random-phase approximation has been made,— quantita- 
tive description of the retardation effect^ still remains a 
serious obstacle for the calculation of T c . 

Recently, the situation has changed by the progress 
in the density-functional theory for superconductors 
(SCDFT).— ~— There, a non-empirical scheme treating 
both the electron-phonon and electron-electron inter- 
actions was formulated referring the Migdal-Eliashbcrg 
(ME) theory of strong-coupling superconductivity.— ~— 
Subsequent applications have shown that this scheme 
can reproduce experimental T c s of various conventional 
phonon-mediated superconductors with error less than a 
few KJJril 

An important advantage of the current SCDFT-based 
scheme is that the particle-hole symmetric component of 
electronic structures around the Fermi level is accurately 
treated on the basis of the density functional calculations. 
In the context of the ME theory, its effect on T c has 
been intensively studied for A15 compounds.— ~— How- 
ever, the particle-hole antisymmetric component ignored 



in the current SCDFT has also been shown to have a sig- 
nificant impact on T c .— Moreover, there has also been 
a great experimental progress of electric-field-induced 
carrier-doping technique, and it has stimulated current 
interest to superconductivity in doped insulators having 
significant particle-hole asymmetry^ - — With this back- 
ground, it is important to extend the applicability of 
SCDFT to systems with particle-hole asymmetric elec- 
tronic structure. 

In Ref. |, the assumption of particle-hole symmetry 
with respect to the Fermi energy (Ep) was originally in- 
troduced to avoid numerical divergence in the exchange- 
correlation kernel. In this paper, we show that we can 
construct an exchange correlation functional for systems 
with asymmetric electronic structure which does not have 
any divergences. In Sec. HTJ we show how the kernel is 
constructed, and discuss how it is numerically stable. 
In Sec. IIII[ we perform calculations of T c for the cases 
where the present improvement becomes indeed impor- 
tant. Section Hvl is devoted to summary and conclusion. 



II. GAP-EQUATION KERNEL FOR 
PARTICLE-HOLE ASYMMETRY 

In this section, we describe the improvement of the 
gap equation based on the current SCDFT to include the 
particle-hole asymmetry effect 
we solve the following gap equation: 
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Here, n and k denote the band index and crystal mo- 
mentum, respectively, A is the gap function, and /3 is 
the inverse temperature. The energy E n \^ is defined as 
£Wk=V Znk + ^nk an d £nk = £nk — M is tne one-electron 
energy measured from the chemical potential fi, where 
e n k is obtained by solving the normal Kolm-Sham equa- 
tion in density functional theory (DFT). 

Functions Z and K. are the exchange-correlation ker- 
nels, which represent the effects of mass renormalization 
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and effective pairing interactions, respectiveiy. In the 
SCDFT, Z is treated by only the phonon contribution as 



because of its antisymmetric dependence on £ nk . The 
resulting form is 
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with Z^k. representing the renormalization of electronic 
states by the phonon exchange, and K. consists of both 
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the electron-phonon /C, ikr[ , k , and electron-electron JC^,^ Further, they modified an unphysically large component 



exchange contributions 
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Formally, Zf^ is derived from the functional derivative 
of the free energy given by the Kohn-Sham perturbation 
theory with respect to the anomalous density^ as 
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Here, 5"k^k' denotes the electron-phonon coupling, and 
/ and I' are defined as 
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with and denoting the Fermi and Bose distribution 
functions, respectively. We note that Zf h originates from 
explicit dependence of the free energy on the anomalous 
density, whereas Z% h comes from the implicit dependence 
via the chemical potential. 

Liiders et al£ reported that this form suffers from nu- 
merical divergence and that the divergent contribution is 
antisymmetric with respect to Ep. Accordingly, in their 
paper they symmetrized Zf h in £ nk , and dropped Z^ 1 



of Z, 
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around E-p, and obtained the following form 
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where J is defined by Eq. (80) in Rcf. |9j. So far, only this 
form has been used in practice. 

In Eqs. © and (fTU)) , the antisymmetric parts of the 

electronic density of states (DOS) and |,9" k !" k k I 2 with re- 
spect to Ep within the phonon energy scale are ignored. 
Within the standard ME theory, the effect of their an- 
tisymmetric parts is assumed to be mino r 20 ' 32 and often 
ignored for simplicity^ 3 - On the other hand, in systems 
with rapidly varying DOS, this treatment is also known 
to lead to a considerable error in the estimate of T c and 
the isotope-effect coefficient^ Below, we will show a nu- 
merically stable form of Z without the symmetrization, 
with which the asymmetry effect is properly treated. 



A. Cancellation of divergent terms 

First, we analytically examine the divergence of 
Z ph,PT . For simplicity, we deal with the nk-averaged 
form of Z ph,PT , defined by 



z^ PT (0=zf(0 + zf(0, 
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where N(£) is the electronic DOS. We also introduce an 
approximation for the electron-phonon coupling 
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where the antisymmetric component of DOS is retained. 
The function a 2 F(u>) denotes the Eliashbcrg function, 
with which the electron-phonon coupling coefficient A 
and the characteristic frequency ui\ n are defined as 
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Starting from the above form for Z, in brief, we show 
that the numerically divergent terms analytically cancel 
with each other so that we can construct a numerically 
stable form. Readers who are not interested in the de- 
tail of this analysis can skip the remaining part of this 
subsection. 

Let us first focus on Zf h . As we show below, the diver- 
gence in this part basically comes from the slow decay of 
order 0(£' -1 ) of the integrands. In order to extract this, 
we transform the following factor appearing in I and I' 
[Eqs. © and ©] as 
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The principal value integral P has been introduced to 
deal with the poles at £'=±w in the right hand side; 
note that the expression in the left hand side docs 
not have these poles. Substituting this expression into 
Eq. (|T2|) . inserting identities 1 = J e?£'<5(£' — £„<k') and 
1 = J dujS(ui — kVk-ic), an d using Eq. (|13j). we obtain 
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Here Ii (i=0, 1) and J, (i=0, 1, 2) denote the terms origi- 
nating from / and J', respectively. From this expression, 
the divergent contribution z^ h ' dlv is extracted as 
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We do not explicitly write the variables (£,£', u;) unless 
necessary in the following. 

Let us examine Zf h ' dlv . For simplicity, we start from 
the case of constant electronic DOS N(^)—N(0) with 
asymmetric energy cutoffs [— L 2 , la] as depicted in Fig.[U 
which is the simplest case of particle-hole asymmetry. We 
then get to 
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Next we carry out the £' integral. Since the integrand is 
of order 0(£' _1 ), integration results in cutoff-dependent 
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FIG. 1: Constant DOS with asymmetric cutoff energies, L\ 
and — L2- 
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Here we have neglected terms of order 0(e _/3w ), 
0(e~P Ll ) and 0(e _/3i2 ), and this approximation is im- 
plicitly used in the following equations with "~" in this 
section. We consequently get to 
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The ^-dependent factor here is reminiscent of Z% h 
[Eq. ©]. In fact, we prove that this term exactly cancels 
-Zf (£) m tne following. 

Next, we consider i?2 h (£) using the approximations of 
the same level. The nk-averaged form is explicitly writ- 
ten as 
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]T(/V2)/cosh 2 [(/3/2)^] ~ 2JV(0). 
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By carrying out the integration J d£ J d£' and ignoring 
terms of order 0[(T/L{) 2 , (T/L 2 )% we get 
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and immediately obtain 
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This expression exactly cancels out Zf h ' div (£) [Eq. ([27)]. 
One can observe the same cancellation for the noncon- 
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FIG. 2: Decomposition procedure of DOS. Typical DOS is 
represented as red-shaded area. 



Here, we approximate the DOS as a step-like func- 
tion with a certain set of energy tics {±£1} (i = 
1, 2, N e ; ei < min[Li, L2]). The approximate function 
is then decomposed into functions 0f(£) = 0(£ — £;) + 
0(-t - a) and 0r(O = 0(£ - £.) - - ei) as 

N® = N^il + ^TlN-e-iO + N+etm}. (34) 
This procedure is schematically depicted in Fig. [5J For 
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, straightforward calculations yield 
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For Z\ h , terms proportional to N^N^ 1 (any double sign) 
seemingly become nonzero, but careful calculations show 
that these terms cancel out each other (See Appendix 
|A")) . Finally, we see that the cutoff-dependent part in Z% h 
exactly cancels Eq. (f3"5"j) . By taking limit for the number 
of tics N € — ^ oo, the present cancellation is proved even 
for general DOSi^ 

Now let us point out the important aspect of Eq. ([Til : 
With the present approximation, the cutoff dependent di- 
vergence cancels if analytically calculated, and the appar- 
ently divergent contribution is in fact nonsingular with 
respect to £, since it is proportional to a nonsingular 

function | - cosh[(/3/2) f ] / s i n hK/3/2)c| ' In ordcr t0 verif y this ' 
we have calculated Z ph ' PT (£) [Eq. (fTTj) ] with a certain 
asymmetric model DOS [depicted later in Sec. IIII] , We 
have generated the energy grid so that it becomes uni- 
form in logarithmic scale. The calculated result is shown 
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FIG. 3: Calculated values of Z ph ' FT {£) [Eq. JIT])] with the ap- 
proximation given in Eq. 1131 Model DOS depicted in Fig. [5] 
with Ef = and model Eliashberg function given by Eq. (|56[) 
with A = 1, we = 400cm" 1 , and a — loe/W were used. The 
calculations were done with T = 0.1K. Logarithmically uni- 
form energy grids with mininum energy cutoff of 10~ 7 (red 
solid line), 10" 8 ( green dashed line), 10 9 (blue dotted line), 
and 10 -10 (violet, dot-dashed line) were used. The DOS val- 
ues for |£|<10T were fixed to iV(0) so that it is consistent 
with the condition (ti^T) assumed in Sec. Ill Al 



in Fig. [3] By lowering the minimum energy cutoff for the 
grid, convergence within order of A has been achieved, 
where the apparent divergence gradually vanishes. How- 
ever, we have also found that the convergence requires 
formidably accurate integration within the energy scale 
smaller than temperature, and so it is extremely difficult 
to achieve in practical calculations. 



B. Z kernel for particle-hole asymmetic systems 

The above analysis specifies the numerically diverging 
(or unstable) but analytically irrelevant terms. In or- 
der to circumvent the numerical difficulty, we propose to 
simply subtract jjP h > dlv and Z\ , obtaining 
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This auxiliary expression shows stable convergence, in- 
cludes the effect of the antisymmetric part of electronic 
states, and agrees with Eq. (jlip with error of order 
0[(T/Li) 2 , {T/L 2 ) 2 \ in the limit where Eq. CTJ]) becomes 
exact. The form given in Eq. (|36p is a generalization of 
the nk-averaged form of Eq. ^ 
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Symmctrization of Eq. (|36|) in £ yields Eq. (|37[) . the lim- 
iting value lim^ [ZP h < aux (£) - ZP h ^ m (0] = because 



of the fact that the integrand for £=0 is symmetric in 
and the two forms show almost the same temperature 
dependence. 

Following the procedure of Ref. [§, we next modify 
ZP h > aux (£) around E F and obtain 
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In practice, we also found that the principal- value in- 
tegral around the singularity at £' = ±u; is numerically 
unstable and its convergence is difficult to achieve with a 
practical computational cost. We therefore introduce an 
even smoothing function p(£' ± u>) whose value is unity 
for \Z'±uj\ > T and of order 0[(£'±w) 2 ] for |f ±u\ < T. 
Substituting p(£' ± ± to) for P[l/(£' ± uj)] and re- 

arranging the terms yield^ 
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We have used p(:r) = [tanh(500/3a;)] 4 in this paper, 
though the calculated value is not sensitive to a specfic 
choice of p(x). 

With the above form, the asymmetry effect is properly 
treated in a numerically stable manner; the smoothness 
of the integrands at f = ±£ ± w,±w, and £ = is re- 
tained. In addition, for approximately symmetric DOS, 
the calculated results are guaranteed to be quite close to 
those with the nk-averaged form of Eq. (jTDJl 
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x[j(z,e,cj)+j(z,-e,Lj)] 

because of the following properties: (i) 
lim 4 _> [ZP h - now (£) - -ZP h (£)] ^ 0, and particu- 
larly for constant DOS, the desirable behavior— 
lim^o ZP h < ncw (£) ~ \im^ Z ph (0 ~ A is satisfied, 
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(ii) temperature dependence of the calculated value is 
similar to that of Z ph (0, and (hi) the symmetric part 
of the integrand, that is, 



2c(€),#c(£) 



J(e,^w)+I(^-^w)-2J(f,^w)-2J-(f,-^w) 



(46) 



agrees with J(£, w)+J(£,— w) in Eq. (|4"5"j) when we im- 
pose a condition uj>\£\~^>T . These properties are demon- 
strated in the following calculations. Further, we have 
found that the nk-resolved counterpart of Eq. (|4T)]) given 
by 
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X [2(Cnk, fn'k', W„k-k') — 2j7(£nk, £n'k', l^Vk-k')] (47) 

is also numerically stable, which is a reasonable general- 
ization of Eq. (fTUj) . Equations (j4T)|) and (|4"T)) are those we 
propose as the new forms of Z for particle-hole asymme- 
try. 



III. ASYMMETRY EFFECT ON T c 

In the rest of the present paper we get insights into the 
asymmetry effect included in the present improvement. 
In Sec. IIII Al we analytically construct a formula of T c 
to discuss in what situation the asymmetry significantly 
affects T c . In Sec. IIII Bl a typical model case is presented 
where our new kernel and the previous one give quite 
different T c , and some remarks on the application to ele- 
mental metals are given. We here note that we focus on 
the nk-averaged form [Eq. (00])], and thereby the present 
scope includes only the asymmetry of DOS, not that of 
the electron-phonon matrix elements. 



A. Analytic calculation 

We first analytically solve the energy-averaged gap 
equation^ 



A(C)=-Z(0A(0 
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with the electron-phonon kernels only. We introduce a 
simple model for which DOS and kernels are given by 
the constant part and the antisymmetric part, given as 

2(0 = MO N(0 = N c (0+N a (0, >C(00) = 

/C c (£,£') + /C a with each part defined within the 
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FIG. 4: (Color online) Schematic description of terms defined 
by Eqs. (jMJ-((ST|. 



Dcbyc frequency wd by 

z c (0=z c , z.(0 = 

N c (0=N c ,N a (0 = 



-Z a , (£<-w') 

o, dei<w') 

-JV., (£<-«') 
0, (fc|<w') 

^a, (£><*/) 



(49) 
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/c c (^,r)=^c, 

[(£ < -w' and £' < -^') 
or (£ > u/ and £' > w'j] 
/C a (£,e'Hir a , [(£<-a/ and £'>o/) 

or (£ > w' and f < -cj')] 

(iei<w' or \e\<uf) 



(51) 



Here, in addition to wd, we have also introduced to' which 
specifies the energy scale where DOS and kernels sub- 
stantially deviate from the values at E-p. These forms 
are depicted in Fig. 0J 

Similarly to McMillan's procedure^ we assume the fol- 
lowing rectangular form for the solution 



'A_, (-a/ > £ > -ud) 
A(0 = {A , (|£|<a/) 

A+, (cj d >£> w') 



(52) 



The condition that a nonzero solution for A exists gives 
T c . By retaining the lowest order term with respect to 
the values with subscript "a" , we obtain 

(l+Z c N a ZM^/u'] Z a 2 ln[o; D ^^ ( ^ 
T c ocexp<^ — at / 1 i ? \ + 7T—^3. f d 53 ) 



iV c (l + Z c ) (l + Z c ) 2 



[K C N C 



Since Z a is dependent on N a through Eq. (|40[) . we can 
substitute an approximate form for Z a in terms of N a /N c , 
and then we obtain 



T c ocexp ■ 



1 + A 



A(r, A)lnr 



(54) 
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Here, r = wd/w', and we have used the following relation 
K C N C = -A and Z c = A (Ref. !). The first term in the 
bracket corresponds to the zeroth order, which appears 
in literatures^ The function A(r, A) is a positive function 
which monotonically increases by increasing r or A and 
converges to a finite value < 1 in the limit r — > oo and 
A —> oo. The detail of the derivation of Eqs. ([5Ti]) and 
([5~4"]) is given in Appendix [B] 

From Eq. (|54|) . we see how the antisymmetric part af- 
fects T c : It reduces T c regardless of the sign of N a and 
its amount becomes substantial when the ratios N a /N c 
or r = lod/lu' are large. In addition, we find the anti- 
symmetric part of the nondiagonal kernel (K a ) does not 
contribute to T c within the lowest order. 



B. Numerical calculation 

The above analysis tells us that the asymmetry effect 
becomes pronounced for the cases where N a /N c and r 
are large. To quantify the effect more explicitly, we cal- 
culated T c by numerically solving Eq. (|48[) using a model 
DOS. We included both the phonon and electron con- 
tribution to the kernels, K. = /C ph + /C o1 and Z = Z ph . 
We employed the nk-averaged form for /C ph [Eq. (23) in 
Ref.llJ, and we treated IC el by an approximate form 



£ el (£,0 



N(0) 



(55) 



with a certain energy cutoff. For Z ph , wc employed the 
previous form [Eq. (|45)) ] and the new one [Eq. (j40|) ] . 

We introduced a step-like DOS as depicted in the up- 
per panel of Fig. [5] so that N a /N c and r can be large. 
This form is characterized by two parameters; the ratio 
of the high-energy and low-energy values N + /N_ , and 
the energy range where DOS changes d. For the Eliash- 
berg function, we used the following model function 



a 2 F(u>) 



A 



2cr, 



zwexp 



uj — ujp, 



(56) 



with width a = cjd/10. We set tup, = 400 [cm -1 ], 
d = Wd/2, and N+/N- = 6. This setting is realistic in 
that one often see a similar dependence in doped semi- 
conductors with quasi-2D Fermi surfaces i 31 ' 35 i 36 

With these settings, we calculated T c for various set- 
tings of Ep. The calculation was performed with a log- 
arithmically uniform energy grid, where the minimum 
energy cutoff was set to |£| = 10 -6 eV and the number 
of points per digit was 20. The energy cutoff for IC cl was 
fixed to 20cV. The parameters A and fi were fixed to 1.0 
and 0.5 for all the setting of Ep so that the resulting T c 
becomes > 10K. 

In the lower panel of Fig. [SJ we plot the ratio of the 
calculated T c using Z ph ' now (£) [Eq. gDJ] to that using 
Z ph (£) [Eq. ([15])] as a function of Ep, where the absolute 
values of T c are given in the inset i2l The ratio is system- 
atically less than unity; the effect of asymmetry lowers 



T c for any Ep, which is consistent with Eq. (|54|) . As the 
Fermi level approaches to the lower edge of the slope, 
the decrease of the ratio becomes significant and finally 
amounts to more than 20%. When the Fermi level is lo- 
cated on the upper edge of the slope, the decrease is not 
appreciable. This Ep dependence can also be consistently 
understood with Eq. (|5"4"|) ; the characteristic energy scale 
(u/) decreases as the Fermi level approaches to the slope, 
and the ratio N a /N c is relatively small when the Fermi 
level is on the upper edge. Thus, the change of T c by the 
present improvement becomes crucial when Ep is close 
to a point where DOS increases rapidly. 

We show the calculated value of Z ph for Ep—— 0.3a>D 
in Fig. |5] to examine the properties of the present form. 
The model DOS is also depicted in the same energy scale 
(green shaded area). For both T=0.01K (thick lines) and 
y^yasym ( thin jjj^ Z P h > ncw (£) becomes asymmetric in 

£, but they are similar to i? ph (£). In particular, the lim- 
iting value for £ — s- agrees with limj^o Z ph (0- These 
features well represent the properties (i)-(iii) described 
in the last paragraph of Sec. Ill Bl We observe a differ- 
ence between Z ph < now (£) and Z ph (£) around |f| > u D , 
which is responsible for the difference of T c asym and T c sym 
in Fig. [SJ Clearly, Z ph ' new (£) becomes larger for the en- 
ergy region with larger density of states, whereas iJ ph (£) 
does not. The dependence of the former is reasonable 
because larger DOS should result in stronger mass renor- 
malization due to stronger total electron-phonon cou- 
pling. Thus we were able to see that the particle-hole 
asymmetric electronic structure is properly treated with 
jjph.ncw ^ n a numer j ca r]y stable manner. 

Finally, let us comment on the application to actual su- 
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FIG. 5: (Color online) (Bottom) Ratio of T c calculated us- 
ing _ZP h ' new [Eq. go])] (T c asym ) to that calculated using Z ph 
[Eq. PS])] (T c sym ) plotted as a function of E F . In the inset, 
absolute values of T c asym and T c sym are shown. (Top) DOS 
used for the calculations. We set N(£)=N(E F ) for |£| > 5eV. 
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FIG. 6: (Color online) Calculated values of Z ph ' nev (£) (red 
solid line) and -Z ph (£) (blue dotted line) for E f =-0.3lu d . The 
green shaded area depicts DOS (arb. unit.) in the same 
energy scale. Note that the values for T = T™ ym (thin lines) 
are multiplied by 0.5. 
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perconductors. For typical superconductors, aluminum 
and niobium, we have carried out stable calculations of T c 
(see appendix [Cj) . confirming again that the listed prop- 
erties (i)-(iii) of Z ph ' ncw are valid. The asymmetry ef- 
fect on T c was, however, estimated to be less than 0.1 % 
for the both systems. This is mainly due to the small 
N a /N c in these systems. The present weak asymmetry 
effect gives support to the validity of the particle-hole 
symmetrizing treatment ^ 13 ' 20 ' 32 for many cases of con- 
ventional superconductors. 



Appendix A: Evaluation of Z ph (£) for nonconstant 
DOS 



We here describe the detail of evaluation of Z| h (£) 
[Eq. ([33]) ] for the nonconstant DOS. Under the de- 
composition approximation of DOS introduced in 
Sec. UTS] [Fig. El Eq. {31])], ZP h < 2 (£) is given by 



IV. SUMMARY AND CONCLUSIONS 



We constructed a new exchange-correlation kernel 
Z ph,ncw in SCDFT to treat the particle-hole asymme- 
try of electronic structure. The obtained nk-averaged 
forms [Eq. (|40[) ] and nk-resolved form [Eq. (|47]) ] do not 
show any divergences or instabilities, and, for the systems 
with good symmetry, well agree with the previous sym- 
metrized forms. By analytically deriving T c formula, we 
found that the asymmetry systematically decreases T c . 
The amount of the decrease becomes substantial when 
the following two ratios are large: The ratio of the anti- 
symmetic component of DOS to the constant component, 
and that of the Debyc frequency to the characteristic en- 
ergy scale of the DOS variation. We also calculated T c 
with a model step-like DOS, showing that the amount 
of the reduction of T c can be more than 20%. With the 
present work, we successfully reinforced the theoretical 
foundation of SCDFT to extend its applicability to sys- 
tems with significant particle-hole asymmetric electronic 
structure. 



Z ph,2 (e) 
1 



p/2 



X cosh[(/3/2)£]sinh[(/3/2)£] 



d£>\i+YyNTej{t')+Npf{0]\i'{U',u)- 

(Al) 



With omitting the temperature-dependent terms of or- 
der 0[(T/e 4 ) 2 ], 0(e-^), 0(e~^ Ll ) and 0(e~' 3L2 ), the 
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integral for each term is carried out as 

Ll dt f Ll dfetftef 

-L-2 J—L 2 

(Li+e l +uj)(L 2 + ej+uj) 



±ln 

Li 



(Li+ej+Lu)(L 2 + e l +uj) 



(A2) 



-L-2 J —L; 

(L 1 +e l +uj)(L 1 + e J +m)(L 2 + e l +Lu){L 2 + e :j +uj) 



±ln 



(L 1 +L 2 +U)) 2 (e i +e j +c J j) 2 



Ll at [ Ll dt'erw&t'^) 



L 2 

In 

Li 

-L 2 

In 

Li 
L-2 

In 

Li 

-L-2 

In 



(Li+L 2 +w) 2 (e. t +w) 2 



{Li+u){Li+ei+uf)(L 2 +u)){L 2 + ei+uj) 



-L 2 



{Li+e i +Lu)(L 2 + e l +Lo) 



(ii+L 2 +w) 2 

d£, f Ll dee+(oi'(t,e,u), 

(L 1 +e l +uj)(L 2 +uj) 



{L 2 + e i +uj){L 1 +uj) 

d£ I 1 dfef (?)!%?,<*>), 



-L 2 



L 2 +€ t +UJ 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 



Using these relations, Z^iS,) is transformed to the same 
form as Eq. (|3"5|) with the opposite sign. 



Appendix B: Derivation of Eq. (|53|) and Eq. (|54|) 



Let us first plug the kernels and the gap function given 
in Eqs. (|49]) -(f52" j) into the energy-averaged gap equation 
[Eq. (HE])]. The condition to have a nonzero solution is 



detM=0, 



Me 



l + Z r 



- Z a + i (N c - N a ) (K c - K a )q N c K c p 
±(N c -N a )K c q l + Z c + N c K c p 



UN C -NJ(K C +K a )q 



N c K cP 



±(N c + N a )(K c +K a )q 
\{N c + N a )K c q 
l + Z c + Z a + \{N c + N a )(K c - K a )q 



(Bl) 



with p and q defined by 

f ' tanh[(&/2^ p , tanh[(^2]g] m<? . 

p= 7o f 7 - f 

Here /3 C denotes the inverse of T c . In order to treat the 
lowest-order contribution, we keep the terms up to the 
second order with respect to the values with subscript 

"a" as 

detM=0 

~(l + Z c f + [K c N c C-K a N c q]{l + Z c f 

+ [-K c K a N?qC-K c N a Z a q-Z 2 a ](l + Z c ) 
-K c N c Zlp, (B3) 



where C = p + q has been introduced. Retaining the 
lowest-order terms, we obtain 



1 + Z N a Z a q Zlq 

K C N C N c (l + Z c ) (1 + Z C ) 2 ' 



Using C = In [2c 1 ojd / (irT c )] with 7 being the Euler con- 
stant, we get Eq. (|53")) . 

Next we consider the relation between Z a and N a . Let 
us start from our newly developed form [Eq. (|40p j. The 
antisymmetric part of the integrand in £' can be written 
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«) - -w) - «) + -w) 

(B5) 

Using the facts I(— —w) = !(£,£', w) and 

i7(— £, — £'? — w ) = 0nc can easily find the 

above part divided by tanh[(/3/2)£] is antisymmetric in 
£. The symmetric part [Eq. (gBJ)] divided by tanh[(/3/2)£] 
is, on the other hand, symmetric in £. Consequently, the 
antisymmetric part of iJP h ' new (£) is yielded by only the 
antisymmetric part of iV(£'). 

Substituting JV(£') = N c (£') +N a (£'), for the antisym- 
metric contribution, we obtain 

= 2 W C J <[(antisym.)], (B6) 

where (antisym.) represents the terms in Eq. (|B5|) . This 
integration is easily performed since we can assume £ 3> 
T, £' > T, and w > T, so that we get 







N cJ 






~ N c 


i 


1 


£ + CJ' + cj £ - 


1- w D + w 



(B7) 



Finally, by using the Einstein spectrum a 2 F(uj) 

i 



■^wd<5(w — wd), we obtain 



|~^ph,ncw ^ph.ncw^ ^ j 



sgn(£)A 



wd(wd - c/) 



(B8) 



where the sign function comes from lirng-^oo 
l/tanh[(/3/2)£]. 

Subsequently, we consider to approximate this form 
as the rectangular form [Eq. (|4"5|) ]. Using the value at 
£ = u>Y) yields 



Z a =X 



1 



N c 3(2r + 1) 



(B9) 



where t=o;d/w'>1. By using this form in Eq. (|53[) and 
substituting K a N c = —A, Z c = A, wc get to 



T c ocexp 




Xh(r) f Xh(r) 
1 + A ~^ 1 + A 



In?' 



(BIO) 



This expression is equivalent to Eq. (|54| , where the func- 
tion A(r, A) is defined as the part in the square bracket 
in Eq. (|B10|) . Its positiveness, monotonic dependences 
and the convergence in the limit r — > oo and A — > oo 
referred in Sec. IIII Al are easily confirmed with this ex- 
pression. We also note that A(r, A) has these properties 
for any values of £ (wd>^>w') in Eq. (|B8|) when we get 
an approximate form of Z a . 



Appendix C: Ab initio calculation 

We applied our formalism to typical weak- and strong- 
coupling superconductors, Al and Nb. We solved the gap 
equation [Eq. ([1}] with the energy-averaged approxima- 
tion for the phonon kernels 



^nk,n ; k' — (£nk}£n'k') 
-£nk=-£(£«k), 



■K2 



'k'> 



(CI) 
(C2) 



where /C ph (£„ k , £n'k') and /C^ kn , k , are defined by 
Eq. (23) in Ref . QJ and Eq. (13) in Ref. HI respectively. 
For the diagonal part Z(£„ k ), we employed -Z ph,now (£nk) 
[Eq. gO}] and Z ph (£„ k ) [Eq. g5)]. We used an accurate 
random-sampling scheme given in Ref. l35l . and the result- 
ing sampling error was not more than a few percent. The 
detailed condition of the calculations is given in Ref. [39|. 

In Fig. the calculated DOS for Al [(a)] and Nb 
[(b)] are shown. The calculated values of Z„ k with the 
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FIG. 7: (Color online) Calculated DOS for Al [(a)] and Nb 
[(b)]. Z kernels for Al [(c)] and Nb [(d)] calculated from 
Eq. (J40J) (solid line) and Eq. (gH) (dashed line): The in- 
put Eliashberg function (dotted line) and DOS (green shaded 
area) are given in the same energy scale, with which A and 
cjin were estimated for Al (Nb) as 0.417 (1.267) and 314 K 
(176 K). 
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two forms are also shown in panel (c) and (d) together 
with the input Eliashberg function a 2 F. For the both 
systems, the DOS (green shaded area) has small but 
nonzero antisymmetric component in the corresponding 
energy scale. However, the two forms of Z yield approx- 
imately the same value. For the calculated T c , regard- 
less of the form of i>„k, we obtained 0.7 and 8.7 K for 



Al and Nb, respectively. Although the gap values calcu- 
lated with iJ ph ' new (£nk) were approximately 0.1% smaller 
than those calculated with -Z ph (£„k) for each set of the 
sampling points, this difference was within the sampling 
error. The present result indicates that the asymmetry 
effect is not significant for T c in Al and Nb. 
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